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NOMENCLATURE 

a  tip  vortex  core  radius 

A  area  of  vortex  sheet 

e(,  angle  of  attack 

c  chord  of  foil 

Cd  drag  coefficient 

Cf  skin  friction  coefficient 

Cp  pressure  coefficient 

Cl  lift  coefficient 

01  induced  drag 

Eln  energy  in  vortex  core 

Eo  energy  surrounding  vortex  core 

Et  total  energy  in  vortex  system 

0,0,  K  circulation 

circulation  at  midspan 
pc  core  circulation 

y  circulation  density  of  vortex  sheet 

F  force  on  cross  section  of  core 

H  ratio  of  displacement  to  momentum  thickness 
in  a  boundary  layer 

T  vector  element  of  vortex 

p  pressure 

r  radial  coordinate  in  (x,r,#)  system 
R  propeller  radius 
Rh  hub  radius 

/ 0  density  of  fluid, dummy  r  coordinate  In  propeller 
calculations 

?  vector  from  control  point  to  vortex  element 
s  span 
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Sn  point  N  on  lifting  line  midway  between  vortices 
t  time 

u  velocity  in  x  direction 

v  velocity  in  y  direction  in  (X/y,z)  system 

velocity  in  r  direction  in  (x/r/0)  direction 

w  velocity  in  z  direction  in  (x,y,z)  system 

velocity  in  direction  in  (x/r,0)  system 

Uo  free  stream  velocity 

ua  x  component  of  disturbance  due  to  a  unit  vortex 
ur  radial  component  of  disturbance  due  to  a  unit  vortex 
ux  tangential  component  of  the  disturbance  due  to  a  unit  vortex 
Uv  velocity  induced  by  vorticity  at  point  x 

V  velocity  vector 

Vo  relative  flow  past  a  blade  section/  neglecting 
disturbance  velocities 

£ 

V  relative  flow  past  a  blade 
JO-  angular  velocity  of  propeller 

Ui  volume  density  of  vorticity  distribution 
7^  dummy  coordinate  for  foil 
|  dummy  x  coordinate  for  propeller 

x  axial  component  of  (X/y/Z)  system  or  (X/T/©)  system, 
y  y  component  of  (X/y/Z)system/  along  span  of  foil  or  blade 
z  other  component  of  (xry/Z)  system 
Ym  y  coordinate  of  midspan 

W  weighting  factor  for  changing  vortex  spacing 
f  velocity  potential 


INTRODUCTION 


The  subject  of  this  report  is  the  interaction 
kinematics  and  dynamics  of  the  trailing  vortex  sheets  in  the 
wake  of  a  propeller.  The  effects  considered  to  be  of  primary 
interest  are  tip  vortex  cavitation/  the  downwash  on  the 
blade  surface/  and  the  resultant  change  in  the  assumed 
loading  and  field  point  velocities  and  pressures  due  to  the 
configuration  of  the  field  of  trailing  vorticity.  It  was 
apparent  early  in  this  work  that  two  important  phenomena 
were  taking  place  in  propeller  wakes  which  were  not 
accounted  for  in  the  classical  propeller  design  models.  The 
trailing  vortex  sheets  were  not  remaining  on  helicoidal 
sheets  but  were  rolling  up  around  a  core  containing  the 
vorticity  shed  from  the  tip  region  and  the  velocity  field  at 
the  propeller  plane  as  a  result  was  not  that  predicted  by 
theory.  In  this  highly  rotational  tip  region/  several 
factors  should  probably  be  taken  into  account  that  are 
normally  considered  negligible  in  potential  flow  theory  as 
applied  to  lifting  surfaces.  There  may  be  strong  streamwise 
pressure  gradients  due  to  either  diffusion  of  vorticity  or 
to  absorbtion  of  the  sheet  vorticity  into  the  core.  These 
considerations  are  especially  important  if  the  phenomenon  of 
tip  vortex  cavitation  is  to  be  understood. 
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Jn  Chapter  I  the  classical  vortex  sheet  rollup  problem 
from  a  wing  with  elliptical  loading  is  considered  in  a 
manner  similar  to  that  of  Westwater  (  1  ).  The  problem  is 
approached  by  a  discrete  vortex  method  modified  by  analytic 
calculations  near  the  ends  of  the  sheet  where  a  square  root 
singularity  exists  in  the  vorticity  distribution,  it  is 
found  that  a  strip  theory  model  is  adequate  for  the  induced 
velocity  calculations  saving  a  great  deal  of  computer  time. 
The  argument  is  made  that  an  infinite  vorticity,  or 
equivalently  a  finite  line  vortex  cannot  exist  in  a  real 
fluid  and  the  tip  vorticity  is  spread  out  over  a  finite 
thickness. 

In  Chapter  II  a  model  of  the  tip  vortex  core  is 
examined  with  the  intention  of  calculating  pressures  in  the 
core  and  the  effect  of  axial  velocity  gradients  In  the  core 
on  the  velocity  field  outside  the  core.  This  analysis  is 
based  on  the  balance  of  pressure  forces  due  to  vorticity 
inside  and  surrounding  the  core  and  the  resultant  inward  and 
axial  velocities. 

Chapter  III  is  an  application  of  the  results  of 
Chapters  I  and  li  to  the  case  of  an  elliptically  loaded 
lifting  line.  The  pressures  in  the  tip  vortex  and  kinematics 
of  the  sheet  motion  are  calculated.  The  velocities  Induced 
at  the  lifting  line  are  calculated,  and  conclusions  drawn  as 
to  how  this  would  alter  the  assumed  load  distribution. 

Chapter  IV  is  an  application  of  the  theory  developed 
earlier  to  a  wing.  The  actual  load,  tip  vortex  pressures. 
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and  vortex  sheet  kinematics  are  calculated  and  compared  with 
experimental  results. 

Chapter  V  concerns  the  application  of  the  vortex  sheet 
kinematics  and  tip  vortex  core  dynamics  to  the  helicoldal 
geometry  of  propellers. 
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CHAPTER  I 

The  instability  of  a  vortex  sheet  shed  from  the 
trailing  edge  of  a  wing  has  been  the  subject  of  much 
contemplation  and  calculation  since  the  beginnings  of 
airfoil  theory.  Nevertheless  for  most  design  applications 
the  trailing  vortex  sheet  is  assumed  to  lie  in  a  flat  plane 
extending  to  infinity  downstream.  Prandtl,  Westwater  and 
Kaden  (^1^3)  among  others  worked  on  the  problem  of  vortex 
sheet  rollup  with  some  success.  The  problem  has  not  received 
much  attention  in  recent  years  because  the  velocities 
induced  on  the  surface  of  a  wing  are  not  greatly  changed 
when  the  sheet  is  considered  to  be  rolled  up  instead  of 
straight.  The  rollup  takes  place  far  enough  behind  the  wing 
to  be  ignored  in  aircraft  applications  in  general.  In  the 
case  of  heavily  loaded  marine  propellers,  however,  the 
effect  may  not  be  negligible.  The  trailing  vortex  sheet  does 
not  proceed  directiy  downstream  but,  as  seen  from  a 
propeller  blade,  rotates  downstream  on  helicoidal  surfaces. 
If  the  tangential  velocities  due  to  rotation  of  the 
propeller  are  large  compared  to  the  forward  velocity  of 
the  ship,  the  sheets  may  pass  very  close  to  the  following 
blades.  The  present  author  has  reexamined  the  wing  problem 
by  several  numerical  techniques  using  a  discrete  vortex 
representation  for  the  distribution  of  vorticity  in  the 
trailing  vortex  sheet. 

The  first  method  chosen  was  to  select  an  elliptical 
loading  on  a  lifting  line  and  divide  the  trailing  vortex 
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sheet  Into  N  discrete  vortices  at  equal  spacing  along  the 
wing.  The  elliptical  loading  Is  chosen  because  this  gives  a 
constant  downwash  along  the  span.  Thus  any  rollup  Is  due  to 
the  Instability  of  the  flat  sheet  model  rather  than  to 
differences  In  velocities  normal  to  the  sheet  along  the 
span. 

The  wing  loading  and  trailing  vorticity  distribution 
appear  below. 


The  velocity  induced  on  the  sheet  by  the  trailing 
vorticity  is  deduced  by  the  law  of  Biot-Savart  (H  ) 


* 
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to  start  out  straight  downstream  from  the  trailing  edge  and, 
as  a  first  estimate  for  this  position,  is  assumed  to  remain 
a  flat  sheet  to  infinity  downstream.  The  sheet  is  divided 
into  discrete  vortices  at  an  even  spacing  on  the  lifting 
line.  (FIG.  1.1  ) 

They  are  also  divided  lengthwise  into  sections  of 
length  delta  X.  The  strength  of  the  trailing  vortex  at 
position  N  is  the  difference  in  bound  circulation  between 
space  position  SI  and  S2  since  a  vortex  may  not  end  in  a 
fluid. 


We  thus  start  the  calculation  with  a  flat  grid  of 
vortices  and  downstream  positions.  The  velocity  induced  by 
the  vorticity  of  the  rest  of  the  sheet  is  then  computed  at 
each  vortex  at  each  discrete  position  downstream  as  follows: 

A(x-  -  X<n  }n*i 

"  /m 

Z*,<|  -  L 

5x  51 

5y  =  /m>  -  ym>l 

~  Z^tP~z  ,TM 

I5I1  *  Sx1  +5i 

The  velocity  at  position  P(n,p)  due  to  a  vortex 
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For  each  position  the  velocity  due  to  each  element  Is 
calculated  and  summed  to  obtain  the  resultant  induced 
velocity  at  each  point  on  the  sheet  due  to  the  rest  of  the 
sheet.  A  matrix  of  induced  velocities  is  formed  due  to  the 
present  assumed  position  of  the  sheet.  The  streamlines 
containing  the  vortex  lines  are  then  found,  starting  at  the 
trailing  edge  with  a  straight  line  cross  section  and 
assuming  that: 

~  ^  +Vn}f>- 1  At 

A  new  sheet  is  thus  formed  by  integrating  the 
streamlines  back.  If  the  flat  plane  wake  were  stable  the 
downwash  would  remain  constant  at  the  trailing  edge  and 
increase  to  twice  this  value  at  cross  sections  far 
downstream.  The  sheet  would  ,  therefore,  remain  straight  In 
cross  section  but  slope  with  decreasing  curvature 
downstream.  Due,  however,  to  truncation  errors  in  the 
numerical  methods  and  to  roundoff  error  in  the  machine,  the 
calculated  downwash  is  never  exactly  constant  across  a  cross 
section.  The  sheet  will  therefore  not  remain  exactly  flat. 
If  it  is  indeed  unstable  in  the  flat  configuration  and 
stable  in  the  rolled  up  position,  as  claimed  by  Kaden,  it 
should  proceed  to  the  rolled  up  position  post  haste.  After 
several  iterations  the  sheet  should  be  rolled  up  into 
whatever  stable  position  it  finds  most  comfortable.  The 
drawing  in  Figured. 2  )  shows  the  results  of  doing  such  a 
calculation  on  a  wing  of  twenty  foot  span  with  an  elliptical 
loading.  For  this  sample  calculation  31  vortices  were  used 
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and  the  calculation  was  carried  far  enough  downstream  to  be 
rolled  up  although  the  spacing  Is  too  wide  to  distinguish 
details  of  the  roll  up  mechanism  In  the  tip  region.  Two 
further  plots  are  shown  in  Figures  (1.3)  and  (1.4).  These 
assume  a  wing  of  two  foot  span  with  an  elliptical  load 
distribution.  The  magnitude  of  the  circulation  at  the 
midspan  is  three  feet  squared  per  second.  The  first  plot  Is 
of  results  using  only  nineteen  trailing  vortices.  It  is 
obvious  that  one  foot  downstream  the  sheet  Is  no  longer 
flat,  but  the  effect  is  far  more  marked  in  the  second  plot 
for  which  99  trailers  were  used.  It  is  obvious  that  the 
choice  of  vortex  spacing  has  an  effect  on  the  ultimate 
shape  computed  in  the  tip  region,  although  the  motion  of  the 
sheet  in  general  away  from  the  tip  may  be  similar  in  the  two 
cases.  In  fact,  the  smaller  the  spacing  is  chosen,  the  more 
dramatic  the  results  at  the  tip.  Since  the  kinematics  and 
dynamics  of  the  fluid  in  the  tip  region  are  of  major 
importance,  especially  if  cavitation  is  to  be  accurately 
anticipated,  this  ambiguity  in  detailed  results  must  be 
resolved.  This  problem,  will  however  be  deferred  for  the 
moment,  to  consider  the  more  mundane  question  of  calculation 
time.  The  two  calculations  for  99  and  19  vortices  consume 
respectively  20  and  12  minutes  of  computer  time  on  a  360 
Model  65.  This  fact  dooms  further  examination  along  this 
path  to  financial  disaster.  No  practical  propeller  design 
method  could  ever  be  sucessful  with  so  expensive  a 
calculation  at  its  base,  considering  the  inevitable  added 
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calculation  required  to  apply  a  method  to  the  helicoidal 
geometry  of  propeller  wakes. 

Method  Two 

A  more  efficient  means  of  calculating  the  trailing 
vortex  paths  was  sought  to  decrease  the  number  of 
computations  required.  Since  the  effect  of  one  element  of  a 
vortex  on  the  velocity  at  a  control  point  is  proportional  to 
the  Inverse  of  its  distance  to  the  control  point  squared/  it 
would  seem  that  the  major  effect  is  that  of  the  nearby 
elements.  A  slight  inaccuracy  in  the  locations  of  the  vortex 
elements  far  from  the  point  in  question  may  not  cause  a 
significant  error  in  the  final  result.  A  strip  theory  model 
was  ,  therefore  ,  used  assuming  that  the  shape  of  the 
sheet  cross  section  at  each  downstream  position  was  constant 
from  that  section  upstream  to  the  foil  and  downstream  to 
infinity.  This  greatly  simplifies  the  calculations  since  the 
trailers  may  be  considered  straight  lines  parallel  to  the  X 
axis  at  each  position  and  the  induced  velocity  may  be 
calculated  analytically. 
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The  procedure  Is  to  start  out  at  the  trailing  edge  with 
the  sheet  straight  in  cross  section  as  before  and  calculate 
the  induced  velocities  on  the  sheet.  The  vortex  locations  at 
the  next  station  downstream  are  then  found  from: 

=*>  +V  At 

where  delta  T  is  the  spacing  in  the  downstream  direction 
divided  by  the  wing  velocity  Uo.  this  procedure  eliminates 
the  need  for  iteration  and  storing  large  arrays  of  positions 
and  induced  velocities  and  gives  results  approaching  those 
of  the  first  method.  See  Figures  (1.5)  and  (1.6) 

The  obvious  fallacy  in  the  methods  discussed  so  far  is 
that  a  discrete  vortex  model  of  the  elliptical  loading  may 
not  give  a  constant  downwash  except  in  the  case  of 
vanishing  spacing.  The  singularity  in  trailing  vortex 
strength  is  represented  rather  poorly  with  such  a  model.  If 
the  model  is  such  that  the  tip  vortex  immediately  moves  out 
of  the  plane  of  the  rest  of  the  sheet,  the  rollup  will 
proceed  immediately.  See  below. 

J  (*° 

If  ,  however  ,  this  area  at  the  end  of  the  wing  is 
represented  in  a  more  mathematically  satisfying  manner  by 
evaluating  the  downwash  at  the  tip  due  to  the  singularity 
analytically,  the  sheet  still  rolls  up.  This  argument  is 
necassary  only  In  an  ideal  fluid  since  in  a  real  fluid  the 
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vortlcity  may  not  remain  infinite  for  any  finite  length  of 
time.  A  line  vortex  may  not  maintain  a  zero  core  thickness 
but  will  diffuse  over  a  finite  area.  Nevertheless  the  ideal 
fluid  argument  is  worth  pursuing  both  for  the  sake  of 
consistancy  and  to  see  if  it  invalidates  the  previous 
results . 

First/  an  expression  is  needed  for  the  downwash  at  the 
tip  due  to  a  loading  in  the  tip  area  of  form:  VfT 


The  velocity  induced  at  a  point  s  near  the  tip  by  such  a 
loading  between  y=0  andy*L  can  be  calculated  from  the  law  of 


Blot-Savart 
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The  point  of  this  discussion  is  that  there  is  a  finite 
downwash  due  to  the  inverse  square  root  singularity  in  the 
trailing  vorticity  distribution  even  when  the  tip  is 
approached.  To  represent  this  effect  with  a  discrete  vortex 
model  it  Is  necessary  to  consider  the  last  vortex  out  to 
experience  a  downwash  due  to  the  singularity  in  the 
vorticity  distribution  in  the  last  interval.  This  problem 
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should  be  taken  account  of  In  the  computation  of  the  motion 
of  the  sheet.  The  method  used  was  to  introduce  a  control 
point  just  inside  the  tip  vortex.  The  strength  of  the  tip 
vortex  is  weighted  so  that  the  proper  velocity  is  calculated 
at  this  point  B. 
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The  velocities  and  positions  of  this  point  B  are  calculated 
at  each  control  point  rather  than  those  of  the  tip  vortex. 
The  tip  vortex  position  is  obtained  by  extrapolating  the 
position  of  the  next  vortex  in  through  point  B. 

If  this  value  Is  used  in  the  induced  velocity 
calculations  according  to  the  strip  theory  of  Method  two, 
the  triggering  is  less  sudden  but  the  final  rollup 
configuration  is  the  same.  See  Figures  (1.7)  and  (1.8) 

Method  Four 

The  problem  of  ambiguous  results  in  the  wing  tip  region 
must  now  be  considered.  Since  the  results  depend  on  the 


CAUCHY  PRINCIPAL  VALUE  AT  THE  TIP 
COMPARISON  WITH  METHOD  1 

FIGURE  1.8 
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vortex  spacing  selected/  there  must  be  a  length 
characteristic  of  the  flow  which  would  give  a  scale  for  the 
spacing  selection.  There  is  no  ambiguity  in  the  results  if  a 
real  foil  is  towed  in  a  real  fluid.  Since  the  problems  occur 
at  the  tip  of  the  wing/  ?t  is  suspected  that  the  length 
scale  sought  has  some  relation  to  the  configuration  of  the 
flow  in  a  real  fluid  in  the  tip  region.  The  most  obvious 
approach  to  this  problem  is  to  define  a  size  for  the  viscous 
core  at  the  tip.  The  usual  means  of  doing  this  is  by 
equating  the  total  kinetic  energy  per  foot  downstream  ?n  the 
ultimate  rolled  up  wake  to  the  induced  drag  on  the  foil. 
Assuming  that  the  final  configuration  is  two  Rankine 
vortices  of  strength  G/the  kinetic  energy  may  be  calculated 
as  follows.  The  vortices  are  assumed  to  be  s  feet  apart  and 
to  have  core  radius  a.  (6) 


The  velocity  potential  of  the  flow  outside  the  core  is  given 
by: 

r .  ± 

is  L  y-  s'  y *y  J 


The  kinetic  energy  outside  the  cores  is  given  by: 
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The  kinetic  energy  inside  the  cores  *s: 


2.TT  t  dr 
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This  gives  a  total  kinetic  energy  per  foot  or  induced  drag 
of : 


since  a  is  presumed  much  smaller  than  the  span: 


if  the  loading  is  elliptical: 


a=  . i“? 7  s 
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This  gives  a  measure  of  the  core  size.  However  there  are  two 
major  objections  to  this  procedure.  First  the  area  close  to 
the  wing  is  of  major  practical  importance  and  the  total 
circulation  certainly  may  not  be  considered  absorbed  into 
the  core  immediately  on  departure  from  the  trailing  edge. 
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Second,  In  a  model  which  assumes  small  pertubation 
velocities  in  the  downstream  direction  any  concentration  of 
vorticity  from  elsewhere  in  the  sheet  in  the  core  Is 
imposs’ble.  The  vortices  must  move  wi th  the  fluid  and  if 
there  is  no  sink  at  the  core  center,  there  can  be  no  net 
flow  into  the  core  region.  By  observing  the  paths  of  the 
vortices  ir,  the  models  considered  so  far.  It  is  obvious  that 
this  is  true.  There  is  no  concentration  of  vorticity  In  the 
core  region  but  rather  a  rollup  of  large  diameter  about  It. 
This  fact  leads  to  the  third  objection  which  Is  that  in  real 
flows  a  concentration  in  the  core  region  is  possible  and  , 
in  fact  ,  likely  by  allowing  axial  changes  in  velocity  and 
pressure  in  the  core.  This  introduces  a  sink  or  source  at 
the  core  center  in  the  two-dimensional  models  and 
invalidates  the  energy  argument. The  method  used  fcr  this 
analysis  must  look  more  closely  at  the  formation  of  the  core 
near  the  wing.  It  was  suggested  by  McCormick  (7  )  that  the 
important  dimension  was  the  boundary  layer  thickness. 

As  was  mentioned  above,  the  deliberations  about  the 
principal  value  of  the  induced  velocity  Integral  at  the  tip 
when  the  slope  of  the  bound  circulation  strength  is  infinite 
are  somewhat  beside  the  point  in  a  real  fluid.  Such  a  strong 
vorticity  as  would  be  expected  at  the  tip  of  a  wing  with 
elliptical  loading  must  be  spread  out  over  a  finite  area  in 
a  real  fluid.  If  this  vorticity  is  being  shed  from  a  wing 
tip  ,  it  must  be  spread  out  over  the  boundary  layer 
thickness  of  the  fluid  on  the  under  (pressure)  surface  of 
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the  wing  when  this  fluid  reaches  the  wing  tip.  if  there  is 
seperation  of  the  flow  around  the  tip#  this  area  would  ,of 
course,  be  greater.  The  tip  vortex  is  ,  therefore  ,  not  a 
line  vortex  but  has  a  finite  core  and  the  velocities  due  to 
it  are  finite  everywhere.  The  core  region  describes  a  solid 
body  rotation  about  its  center. 

The  infinite  value  of  vorticity  at  the  tip  has  now 
disappeared  ,  and  the  core  vorticity  is  symmetrical  about 
the  center  of  the  core.  The  velocity  induced  at  the  center 
of  the  core  by  the  core  itself  is  therefore  zero.  The  rest 
of  the  sheet  induces  an  upward  velocity  at  the  core  center. 
The  rollup  phenomena  is  now  easily  visualized.  The  end  of 
the  sheet  starts  up  while  the  rest  of  the  sheet  bends  down. 
The  core  then  pulls  the  sheet  out  and  the  sheet  pulls  the 


represents  physical  reality  best  ,  it  will  be  used  in  the 
remainder  of  the  paper.  Results  using  this  method  are 
similar  to  those  of  the  other  systems  (See  Figure  1.9  ),  but 
for  further  investigations  into  the  dynamics  of  the  tip 


Cross  section  token  at  one  foot  behind  a  lifting  line  with  an  elliptical 
load  distribution.  The  circulation  at  midspan  is  3ft2/sec. 

19  trailers  are  used  the  last  3  of  which  are  in  a  boundary  layer  of 
thickness  .05  feet.  The  last  5  trailers  are  at  the  fine  spacing  of  .01667ft 
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vortex  core  region  this  model  seemed  most  appropriate. 

The  vortex  spacing  at  the  tip  is  chosen  small  enough  so 
that  several  vortices  are  included  in  the  tip.  The  strength 
of  the  tip  vortex  Is  then  evenly  divided  between  these 
vortices  that  lie  within  a  boundary  layer  thickness  of  the 
tip.  The  rollup  kinematics  show  little  dependence  on  the 
choice  of  boundary  layer  thickness.  This  question  does 
become  more  important  when  the  dynamics  of  the  core  region 
are  being  considered.  The  minimum  pressure  at  the  center  of 
the  rotating  region  is  for  instance  dependent  on  the  core 
diameter  selected.  The  success  of  this  model  in  eliminating 
the  previous  difficulties  with  ambiguous  results  due  to 


l 

1 

! 


various  choices  of  spacing  may  be  observed  by  comparing 
figures  (1.9)  and  (1.10). 


PAGE  29 


CHAPTER  I ! 

The  preceding  analysis  ,  as  may  be  seen  from  the 
results  of  Chapter  1,  does  not  lead  to  a  concentration  of 
vorticity  in  the  core  region,  although  it  does  give  a  rollup 
of  the  vortex  sheet.  As  a  matter  of  fact  ,  it  is  impossible 
for  any  concentration  of  vorticity  to  occur  in  a  model  which 
assumes  the  streamwise  velocity  to  be  constant.  The  only 
mechanism  which  can  produce  a  concentration  of  vorticity, 
which  implies  a  concentration  of  fluid  since  the  vortices 
move  with  the  fluid,  is  a  sink  at  the  core  in  the  two 
dimensional  strip  theory  model,  which  corresponds  to  an 
increase  in  the  axial  velocity  in  the  three  dimensional 
problem.  Batchelor  in  (  8  )  investigated  the  diffusion  of  a 
line  vortex  and  the  axial  velocities  and  core  pressure 
gradients  associated  therewith.  His  analysis  is  ,  however  , 
limited  to  the  region  far  downstream  from  the  wing  where  all 
vorticity  is  considered  concentrated  in  the  core  region  and 
there  is  no  additional  vorticity  entering  the  core  from  the 
sheet.  His  model  is  ,  therefore  ,  one  with  a  deacceleration 
of  the  axial  velocities  due  to  an  increase  of  pressure 
downstream.  The  pressure  rises  downstream  due  to  the  viscous 
deacceleration  of  tangential  velocities  in  the  core  region 
and  the  resultant  decrease  in  centrifugal  forces.  The  core 
therefore  expands  and  weakens  downstream  in  the  far 
downstream  region.  Professor  Sheila  Widnall  and  Thomas 
McMahon  conducted  a  preliminary  investigation  of  the  tip 
vortex  rollup  problem  using  Batchelor's  model.  (^) 
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Batchelor  begins  his  analysis  with  the  assumption  that 
the  motion  is  axisymmetric  and  the  following  equations  of 
motion  are  valid  in  axisymmetric  incompressible  fiow: 

-I-  U.  +  v  Jl/3r  =  a 

-*■  a  -">r  -  ->  +  »  />V- 
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where: 

U  »J-_  Ji.  +  +  -L  -i- 


The  assumption  is  made  that  axial  gradients  are  small 
compared  to  radial  gradients.  Thus: 

i  <<  sL 

c)x  dr 

V  «  u. 

These  assumptions  are  valid  for  the  tip  vortex  dynamics 
considered  in  this  report.  The  radial  velocity  is  never  more 
than  two  percent  of  the  axial  velocity  and  axial  changes  in 


veloci ty 


small  compared  to  radial  changes.  Thus 


equations  -1-  through  -4-  become: 
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The  integral  form  of  equation  '6-  is: 

-1-  p-- 

r 

In  the  region  far  downstream,  Batchelor  finds  a  solution  for 
this  system  of  equations  of  form: 


If  the  core  radius  a  is  defined  as  the  radius  where  the 
tangential  velocity  w  is  a  maximum,  this  core  radius  may  be 
solved  for: 


This  assumes  that  the  core  had  zero  radius  (was  a  point 
vortex)  at  x=0.  The  change  of  core  radius  with  distance 
downstream  is: 

d  a.  _ 

'lz_  ax  ix 

The  magnitude  of  this  quantity  is  important  in  determining 
the  importance  of  viscous  effects  in  the  core  region  near 
the  wing.  In  the  problems  dealt  with  here  .02  feet  may  be 
taken  as  a  typical  value  for  the  core  radius  near  the  wing 
as  derived  from  the  boundary  layer  thickness.  Uo  is 
typically  25  feet  per  second  and  the  kinematic  viscosity  of 
water  is  about  ten  to  the  minus  five  feet  squared  per 


second.  These  values  give  for  the  downstream  distance  x 
which  is  a  measure  of  how  long  it  takes  a  point  vortex  to 
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diffuse  to  a  core  radius  a: 


H-  v  15X0.OOOH  _ 

}i  H7io=rrm- 


n  8  ft. 


In  other  words,  the  vortex  would  have  to  start  as  a  point 
vortex  198  feet  upstream  to  have  diffused  to  a  core  radius 
of  .02  feet  at  the  trailing  edge.  At  this  value  of  x: 


-14- 
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Therefore,  it  is  concluded  that  in  water  the  core  radius 
does  not  change  fast  enough  due  to  viscosity  to  have  any 
significant  effect  on  results  in  the  region  close  to  the 
wi  ng. 


Near  the  wing,  however,  the  model  is  considerably 
altered.  It  is  possible  to  have  vorticity  from  the  sheet 
enter  the  core,  thus  increasing  the  circulation  and  lowering 
the  core  pressure  downstream.  This  gives  rise  to  an 
increasing  axial  velocity  in  the  core  region  and  an  inward 
component  of  velocity  at  the  core  surface. 

The  mechanics  of  the  analysis  proceed  as  follows.  The 
wing  is  divided  up  into  trailers  as  before  and  the  last  few 
are  considered  to  be  in  the  boundary  layer  and  of  constant 
strength.  The  boundary  layer  is  assumed  to  be  identical  with 
the  core  thickness  and  the  vorticity  is  assumed  to  be 
diffused  over  the  surface  of  the  core.  The  model  is 
therefore,  at  each  downstream  section,  a  Rankine  vortex  with 
a  group  of  trailing  vorticies  outside  it. 
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Initially  at  the  wing  when  the  sheet  is  straight/  the 
minimum  pressure  in  the  core  is  easily  calculated  assuming 
axial  symmetry  and  a  small  effect  of  the  sheet  on  the  core. 
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The  pressure  force  on  the  upstream  end  of  the  core  is: 
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Further  downstream  two  effects  act  to  decrease  the  core 
pressures  and  forces.  The  rollup  of  the  vortex  sheet 
increases  the  circulation  surrounding  the  core  and  any  flow 
Into  the  core  brings  In  circulation  from  the  sheet.  The 
first  effect/  that  of  the  sheet  rollup,  is  small  due  to  the 
fact  that  the  sheet  vorticity  is  small  compared  to  that  in 
the  core  and  the  1/R5  term  in  the  pressure  1  ntegral  (oJoov*) . 
This  effect  is  approximated,  however  ,  in  the  following  way 
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both  for  the  sake  of  completeness  and  to  trigger  the 
concentration  phenomenon. 

The  analysis  assumes  axial  symmetry  and  ,  therefore  / 
the  outside  vorticity  must  be  considered  to  be  evenly 
distributed  about  the  axis.  This  is  approximated  by 
calculating  the  angle  subtended  by  each  section  of  the  sheet 
on  the  core  and  considering  the  effective  strength  of  that 
vortex  on  the  core  pressure  is  this  angle  divided  by  two  pi. 
times  its  real  strength.  Its  radius  from  the  core  center  R 
is  considered  to  be  that  of  the  discrete  vortex. 


This  model  is  ,  of  course  ,  crude  but  the  effect  is  in 
any  case  small  compared  to  the  concentration  effects  in  the 
core  region  itself. 

The  major  effect  is  that  of  the  increasing  core 
circulation  due  to  absorbtion  of  the  trailing  vortex  sheet 
into  the  core.  This  may  be  represented  as  follows.  A  length 
delta  X  of  the  core  is  considered  with  a  pressure  force  F 
acting  on  it  in  the  axial  direction.  A  radial  velocity  V  is 
assumed  at  first  to  be  zero.  Using  the  momentum  equation  in 
the  x  direction  over  a  control  volume  of  radius  a  and  length 
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delta  x: 


X  momenta m  oat  /sc c.  - 

X  momentum  in  /sec.  “^OTTQ-^U.*  +  2.T^<X./,  v  UQ  aa 
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The  original  F  is  calculated  from  the  small  rollup  effect 
described  in  the  previous  pages.  The  resulting  radial 
velocity  V  means  that  a  sink  of  strength  2  aV  per  unit 
length  is  pulling  fluid,  and  sheet  vorticity  with  it  in 
toward  the  core.  The  resulting  Increase  in  core  strength  •» 
calculated  by  assuming  that  all  sheet  vorticity  within  a 
radius  difference  of  V  delta  T  is  added  to  the  core  strength 
and  subtracted  from  the  sheet  strength.  (The  stretching 
effect  is  neglected. 
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This  gives  a  decreased  pressure  on  the  downstream  end 
and  a  resultant  increase  of  axial  velocity  and  sink 
strength.  The  calculation  is  repeated  until  the  radial 
velocity  converges  to  a  limit.  The  new  vortex  positions  may 
now  be  calculated  and  the  process  moved  to  the  next  down 
stream  position. 


it  is  entirely  possible  that  the  tangential  velocity  in 
the  core  is  not  proportional  to  the  distance  from  the  core 
center.  However,  it  may  be  assumed  that  the  tangential 
velocity  is  a  smooth  function  of  radius  from  core  center  to 
core  radius,  if  ,  for  example  ,  a  parabolic  velocity  profile 
is  assumed,  the  pressure  at  core  center  is  only  changed  5.3 
percent.  The  Rankine  vortex  is  therefore  assumed  to  be  a 
reasonable  approximation  for  the  purpose  of  calculating 
minimum  pressures.  It  is  also  assumed  that  the  axial 
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velocity  is  uniform  across  the  core  crosssection  and  equal 
to  Uo  outside  the  core.  The  more  rigorous  method  for 
selecting  these  initial  axial  velocities  would  be  through 
calculating  the  head  at  each  radius  and  completing  the 
energy  balance.  However,  this  fluid  in  the  core  has  passed 
through  the  boundary  layer  on  the  wing  and  has  iost  an 
indeterminant  amount  of  energy.  For  want  of  a  better 
apprcmimation  the  axial  velocity  is  ,  therefore  ,  assumed 
uniform. 
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CHAPTER  l I  I 

The  purpose  here  is  to  apply  the  theory  of  the  previous 
chapter  to  a  lifting  line  load  and  try  to  use  it  to 
calculate  the  resultant  core  pressures  and  the  motion  of  the 
sheet.  It  is  apparent  that  the  minumum  pressure  coefficient 
is  dependent  on  core  diameter,  load  and  the  load 
distribution  over  the  span,  and  the  amount  of  rollup  that 
has  taken  place. 

To  evaluate  these  effects,  a  lifting  line  with  a  span  of 
two  feet  was  loaded  with  various  elliptical  loads  and 
boundary  layer  thicknesses.  The  motions  of  the  trailing 
vortex  sheets  are  calculated  and  the  pressures  in  the  core 
region  estimated.  The  greatest  attention  was  given  to  the 
wake  near  the  trailing  edge  of  the  foil  since  this  has  the 
most  effect  on  the  foil  itself. 

The  program  operates  in  a  manner  similar  to  that  of 
method  four  in  Chapter  I.  The  only  modification  necessary  is 
a  line  sink  at  the  core  center  and  axial  and  radial 
velocities  are  calculated  in  the  core  region.  The  minimum 
pressure  is  calculated  at  the  core  center  at  each  down 
stream  position.  A  flow  diagram  for  this  program  is  shown  in 
Figure  (3.1  ).  Runs  were  made  at  a  variety  of  different 
boundary  layer  thickness,  maximum  circulations,  and 
downstream  spacings  to  determine  the  effect  of  these 
parameters  on  the  minimum  pressure. 

In  Figure  (3.2  )  the  minimum  pressure,  y  coordinate  of 
the  core  center,  axial  velocity  in  the  core,  and  radial 


FIG.3.J  LOGIC  FLOW  GRAPH  FOR  VORTEX  SHEET  AND  TIP  VORTEX 
CALCULATIONS  WITH  A  SINK  AT  THE  CORE  CENTER.  A 
STRAIGHT  LIFTING  LINE  WITH  ELLIPTICAL  LOAD  DISTRIBUTION 
IS  ASSUMED. 
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velocity  at  the  core  boundary  are  plotted  for  various 
distances  downstream  for  a  maximum  circulation  of  3  feet 
squared  per  second  and  a  boundary  layer  thickness  of  .01 
feet.  The  results  show  that  there  Is  an  Immediate  radial 
flow  into  the  core  close  behind  the  trailing  edge  with  a 
corresponding  drop  in  the  pressure.  This  is  why  tip  vortices 
are  sometimes  observed  to  cavitate  a  short  distance  behind 
the  foil  rather  than  exactly  at  the  foil.  It  is  apparent 
that  for  the  minimum  pressure  calculations  the  spacing  in 
the  downstream  direction  must  be  taken  quite  fine  near  the 
foil  to  obtain  the  details  near  the  foil.  If  only  the 
kinematics  of  the  sheet  are  required,  the  spacing  may  be 
much  coarser.  In  Figure  (3.3  )  the  minimum  pressure  is 
plotted  versus  maximum  circulation  for  a  boundary  layer 
thickness  of  0.001  feet,  a  thickness  that  corresponds  to  a 
laminar  flow  on  the  pressure  side  of  the  rudder  model  of 
Chapter  IV.  The  minimum  pressure  is  strongly  dependent  on 
the  maximum  circulation  especially  at  very  small  boundary 
layer  thicknesses  as  would  seem  reasonable  according  to 
equation  -9-  of  Chapter  II.  The  vertical  component  of 
induced  velocity  was  plotted  for  a  maximum  circulation  of 
1.86  feet  squared  per  second  on  a  lifting  line  of  1.67  feet 
span.  See  Fig.  (3.4).  It  is  clear  that  the  downwash  from 
this  trailing  vortex  sheet  is  anything  but  constant  although 
the  load  is  elliptical.  The  induced  velocity  is  in  fact 
positive  at  the  tip.  This  clarifies  some  d i screpenci es 
between  propeller  design  theory  and  practice,  assuming  that 
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the  results  are  analogous  for  propeller  blades.  A  blade 
designed  for  elliptical  load  distribution  would  in  fact  be 
much  more  highly  loaded  in  the  tip  region  than  anticipated 
due  to  this  upwash.  The  tip  vortex  would  be  correspondingly 
stronger  and  the  minimum  pressure  in  the  core  would  be 
lower.  Propeller  designers  have  tended  to  reduce  their 
design  load  drasticly  in  the  tip  region  to  prevent  tip 
vortex  cavitation.  These  results  indicate  that  the  load  is 
inevitably  greater  than  predicted  in  the  tip  region  if 
lightly  or  moderately  loaded  propeller  theory  is  used.  The 
boundary  thickness  must  be  known  accurately  if  the  minimum 
pressure  is  to  be  predicted  with  any  confidence.  The 
dependence  of  minimum  pressure  on  boundary  layer  thickness 
is  indicated  in  Figure  (3.5) 


PAGE  44A 


Boundary  Layer  Thickness  in  feet 


MINIMUM  PRESSURE  COEFFICIENT 
AS  A  FUNCTION  OF  BOUNDARY 
LAYER  THICKNESS  AT  TRvAiln& 

EDGE  OF  TIP,  PRESSURE  S  ID  £ 
FIGURE  3.5 
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CHAPTER  IV 

The  results  of  Chapter  III  offer  the  possibility  of 
experimental  verification.  The  minurnum  pressure  coefficient 
may  be  determined  by  observing  the  point  where  cavitation 
begins  in  the  tip  vortex  core  in  a  moving  stream  of  water  of 
known  pressure  and  velocity.  The  calculations  of  the 
preceding  chapter  provide  a  method  for  determining  the 
minimum  pressure  coefficient  in  the  vortex  core.  When  this 
pressure  coefficient  becomes  equal  to  the  cavitation  index 
of  the  flow,  a  vapor  core  should  be  observed  in  the  tip 
vortex.  The  foil  selected  for  tnese  tests  was  a  bronze 
rudder  model  with  the  following  characteristics: 


span 

root  chord 
tip  chord 

sweep  of  1/4  chord 
geometric  aspect  ratio 
taper  ratio 
section  type 

thickness  to  chord  ratio 


10" 

8.925" 

5.360" 

11  degrees  aft 
1.40 
.60 

NACA  66 

root  .2,  tip  .1 


The  boundary  layer  thickness  is  of  course  a  critical 
parameter  in  the  calculation  of  minimum  pressures  in  the 
vortex  core  as  shown  in  Fig.  (3.5).  To  make  an  accurate 
estimate  of  boundary  layer  thicknesses,  it  is  necessary  to 
know  the  pressure  distribution  on  the  pressure  side  of  the 
foil.  This  is  calculated  for  various  angles  of  attack  using 


IV 
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the  published  character i st ics  for  this  two  dimensional 
section  from  Abbott  and  Doenhoff  (  ) .  Of  course  the  interest 
here  is  only  in  the  flow  aft  of  the  stagnation  point.  The 
high  pressures  forward  of  the  stagnation  point  may  be 
disregarded  since  this  fluid  passes  to  the  other  side  of  the 
wing.  For  the  integration  of  the  boundary  layer  equations 
back  from  the  stagnation  point  to  obtain  the  displacement 
thickness  t  am  indebted  to  Professor  Jerome  Milgram  (II)  who 
has  programmed  the  methods  of  TruckenbrodtOW  and  Spence((3). 
In  both  of  these  methods  Ludv/eig  and  Tillmann's  (J3)  equation 
for  the  wall  shear  stress  is  used: 


_  -.2.68 

.146  R&  10 


678  H 


It  Is  assumed  that  the  boundary  layer  is  turbulent  in  these 
calculations.  This  assumption  is  reasonable  for  ship 
propellers  but  leads  to  difficulty  in  the  experimental  work 
with  this  rudder  model.  It  was  necessary  to  add  strips  of 
plastic  aluminum  to  the  pressure  side  of  the  foil  to  produce 
a  turbulent  condition  on  the  model.  Cross  flow  effects  were 
not  taken  into  account  and  this  may  cause  the  calculated 
boundary  layer  thickness  to  be  greater  than  the  true 


thickness.  Sample  result:-* 

are  shown 

in  Figure 

(4.1). 

Truckenbrodt 's  method  does 

not 

converge 

properly 

i  n  an 

adverse  pressure  gradient 

and 

therefore 

only  the 

me thcd 

developed  by  Spence  is  used 

in 

the  low 

angle  of 

attack 

range,  it  is  felt  by  Professor  Milgram  that  the  difficulty 
in  the  application  of  Truckenbrodt ‘s  method  lies  in  his 


Ang/e  of  aJtcw/c  ia  d legrees 
FI  CURE.  4. I 
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numerical  handling  of  the  equations  at  high  values  of  H  (the 
ratio  of  displacement  thickness  to  momentum  thickness) 
rather  than  an  inherent  fault  of  the  method,  in  the  range  of 
values  of  angle  of  attack  above  5  degrees/  the  two  methods 
agree  nicely.  It  is  felt  that  these  methods  give  a  better 
estimate  of  the  boundary  layer  thickness  than  that  used  by 
McCormick  (7),  which  does  not  take  the  pressure  gradients 
due  to  the  potential  flow  around  the  foil  into  account. 

The  computer  program  described  in  Chapter  IV  Is  used  to 
calculate  the  motions  of  the  sheet  and  minimum  pressure 
coefficient  in  the  vortex  cores.  The  loading  is  assumed  to 
be  elliptical  so  that  the  downwash  is  constant  across  the 
span.  For  a  span  of  1.67  feet  (  twice  the  real  span  due  to 
the  image  effect  of  the  tunnel  wall  )  and  a  chord  at 
midsection  of  .745  feet/  the  circulation  at  midsection  is 
computed  as  follows  for  angle  of  attack*.: 


c 


L'N“ 


^>uQr 

j^°U02C 


i  +zc/s 


3.32< 


With  a  tunnel  water  velocity  of  25 


feet  per  second: 


r-.5i7  < 


where °Cis  in  degrees.  The  downwash  may  also  be  computed 
assuming  that  the  sheet  remains  flat: 


W  = 


C|_  c 

tTs~ 
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The  lift  coefficient  for  this  low  aspect  foil  was  checked 
using  a  lifting  surface  computer  program  developed  by  S. 
Widnall(l5).  This  program  gives  a  slightly  lower  lift 
coefficient  at  each  angle  of  attack. 

Cu=  2.5  7 

The  results  both  from  the  program  and  experimental  are 
plotted  in  figure  (4.6).  The  points  represented  by  small 
circles  on  this  plot  were  measured  by  observing  the 
pressure  at  which  cavitation  first  appeared  in  the  tip 
vortex  core.  The  pressure  in  the  cavity  was  then  assumed  to 
be  vapor  pressure.  Some  error  may  be  expected  due  to  this 
assumption.  The  cavity  pressure  may  well  be  higher  due  to 
air  content  in  the  water.  The  consistency  of  the  results  as 
the  tunnel  remained  under  vacuum  for  several  hours,  however, 
indicates  that  this  effect  was  not  major.  The  initial  runs 
were  made  with  a  tip  rounded  with  plastic  putty.  The 
initial  job  was  rough  and  the  maximum  span  was  about  at  the 
three  quarter  chord  point  rather  than  at  the  trailing  edge. 
This  caused  the  tip  vortex  to  form  forward  of  the  leading 
edge  and  to  seperate  around  the  tip.  The  Reynolds  number  of 
the  flow  around  this  model  was  about  10^  and  some  difficulty 
was  encountered  in  obtaining  turbulent  flow  in  the  favorable 
pressure  gradient  on  the  pressure  side  of  the  foil.  The  tip 
was  rebuilt  with  putty  and  the  foil  was  carefully  smoothed 
with  steel  wool.  See  Fig.  (k.5).  Turbulence  stimulation  was 
provided  by  a  spanwise  strip  of  putty  along  one  side  of  the 
foil  at  the  point  of  maximum  thickness.  The  results  obtained 


PAGE  51 


with  this  foil  are  also  plotted  in  Figure  (4.6).  It  is 
evident  that  when  the  smooth  side  is  the  pressure  side  the 
foil  cavitates  much  more  readily  at  low  angles  of  attack 
than  when  the  turbulence  stimulation  is  on  the  pressure 
side.  Transition  to  turbulent  flow  is  indicated  by  the 
results  for  the  smooth  pressure  side  at  high  angles  of 
attack.  The  critical  cavitation  numbers  for  the  two  cases 
then  approach  each  other.  These  results  compare  favorably 
with  those  of  McCormick(7)  although  the  hysterisis  effect 
that  he  observed  was  not  apparent  in  these  tests.  The 
critical  cavitation  number  was  the  same  when  the  critical 
point  was  approached  from  low  pressures  or  from  high 
pressures.  The  results  indicate  that  the  foil  cavitates  at 
slightly  higher  pressures  than  anticipated  by  theory.  This 
may  be  explained  by  the  fact  that  the  lift  in  the  tip  region 
is  higher  than  predicted  by  theory  as  discussed  in  Chapter 
III.  The  laminar  flow  boundary  layer  thickness  is  calculated 
without  taking  account  of  the  pressure  gradient  along  the 
blade.  Since  this  report  is  concerned  with  marine 
propel lers,  which  operate  in  a  highly  turbulent  flow 
normally,  the  laminar  results  are  interesting  only  In  a 
purely  academic  sense.  If,  however,  the  reader  is  tempted  to 
apply  these  results  to  cases  where  the  flow  is  truly  laminar 
it  would  be  advisable  to  program  a  laminar  boundry  layer 
thickness  calculation  for  a  variable  pressure  gradient  in 
the  outer  flow. 
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The  trajectory  of  the  tip  vortex  was  plotted  (Fig. 
4.2)  for  the  rudder  model  at  an  eight  degree  angle  of 
attack.  Its  position  at  one  foot  downstream  from  the 
trailing  edge  is  .135  feet  above  the  center  of  the  sheet. 
This  agrees  with  the  photograph  (Fig. 4.3)  of  the  rudder 
operating  in  these  conditions  in  the  Propeller  Tunnel.  (The 
angle  of  attack  is  negative  in  the  photograph  so  that  lift 
is  down  and  the  tip  vortex  is  the  lower  curve.)  It  was 
observed  during  these  tests  that  the  tip  vortex  moved  in  the 
free  stream  direction/  in  other  words  along  the  axis  of  the 
tunnel/  shortly  after  leaving  the  trailing  edge  while  the 
rest  of  the  trailing  vortex  sheet  curved  downward  as 
predicted  by  airfoil  theory.  This  agrees  with  the 
calculations  of  this  paper  (see  Fig.  4.2)  and  is  explained 
by  the  fact  that/  if  most  of  the  trailing  vorticity  is 
concentrated  in  the  core/  there  is  nothing  to  induce  a 
pertubation  velocity  on  the  tip  vortex  but  its  image,  which 
is  far  away.  The  downwash  due  to  the  image  vortex,  assuming 
that  it  is  1.67  feet  away  and  that  the  circulation  at 
midspan  is  3  feet  squared  per  second  is  less  than  one 
degree.  Therefore  the  tip  vortex  moves  in  the  freestream 
direction  as  soon  as  it  has  absorbed  most  of  the  rest  of  the 
vorticity.  The  spanwise  trajectory  of  the  tip  vortex  Is 
shown  in  the  photograph  of  Figure  (4.5).  The  eight  degree 
angle  of  attack  with  25  feet  per  second  water  speed 
condition  was  maintained  for  this  photograph.  Again  the  tip 
vortex  curves  sharply  in  the  vicinity  of  the  foil  and  then 
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proceeds  nearly  straight  downstream.  The  calculated  spanwise 

;  4 

trajectory  is  shown  for  comparison  in  Figure  C4.7). 


5 


•  Tip  vortex 


FIG.  4.2  calculated  trajectory  of  the  tip  vortex  core 


and  the  wake  at  midspan  for  the  rudder  model 
at  8  degrees  angle  of  attack  and  25  ft. /sec. 
water  speed. 


1  foot 

FIG.  4.3  A  photograph  of  the  rudder  model  operating 

at  8  degrees  angle  of  attack  and  25  ft. /sec. 
water  speed.  The  lower  streak  is  the  tip 
vortex.  The  upper  streak  is  the  rest  of 
the  wake. 


FIG.  4.4  A  top  view  of  the  rudder  model  showing 
turbulence  stimulating  strips. 
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4.5 


A  top  view  of  the  rudder  model  opera 
at  8  degrees  angle  of  attack  and  25 
water  speed. 


ting 

ft. /sec. 
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Spanwise  calculated  trajectory  of  the  tip  vortex  core  for  the  rudder 
model  at  8  degrees  angle  of  attack  and  25  ft. /sec.  water  speed.  The 
ordinate  is  distance  from  the  tip.  This  calculation  is  for  compariso 
to  the  photograph  in  Fig.  4.5. 
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Feet  in  from  tip 
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CHAPTER  V 


If  the  wake  conf iguation,  induced  velocities/  and 
minimum  pressures  due  to  a  lifting  line  wing  can  be 
calculated  as  in  Chapter  IV/  there  is  no  reason  why  the 
lifting  line  model  of  a  propeller  can  not  be  handled  in  the 
same  manner.  The  propeller  is  more  difficult  only  in  the 
sense  that  its  geometry  is  cylinderical  rather  than  flat  and 
the  discrete  trailing  vortex  filiments  will  follow  more  or 
less  helical  paths  rather  than  straight  lines  downstream  as 
seen  by  an  observer  on  a  blade.  If  the  trailing  vortices  had 
no  effect  on  themselves,  the  axial  velocity  would  be  the 
forward  velocity  of  the  ship  or  the  tunnel  water  speed  in  a 
propeller  testing  tunnel.  The  wake  of  the  ship  is  neglected 
although  modifications  to  include  this  are  in  the  planning 
stage.  The  tangential  velocity  would  beilr  where  jTL.  is  the 
rotational  speed  of  the  propeller  and  r  is  the  distance  out 
from  the  shaft  center.  These  axial  and  tangential  velocities 


are  augmented 

by  the  velocities  i 

nduced 

by 

the 

trai 

1  i  ng 

vortices  from 

the  blade  under  con 
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and 
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those  shed  from  all  other  blades. 

1 1 

i  s 
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possible  that  in  the  case  of  a  heavily  loaded  propeller  the 
radial  induced  velocities  are  important,  especially  in  the 
determination  of  the  roll  up  kinematics.  The  vortices  shed 
from  each  of  K  blades,  which  are  assumed  to  be  similar,  are 
assumed  to  travel  in  similar  trajectories,  so  that  only  the 
position  of  one  vortex  sheet  must  be  calculated. 
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The  strip  theory  model  is  employed  as  in  the  wing  model 
at  downstream  stations.  The  r  and  or  equivalently  Y  and  Z 
coordinates  of  the  intersection  of  each  trailer  with  the 
plane  perpendicular  to  the  X  axis  at  x  is  calculated  from 
positions  at  x- ax  as  before. 

The  only  difference  between  this  computation  and  that 
for  the  wing  is  that  each  trailing  vortex  is  now  assumed  to 
travel  in  a  helical  trajectory  rather  than  in  a  straight 
line  downstream.  The  induced  velocity  at  the  vortices  is 
only  weakly  dependent  on  the  pitch  angle  of  the  helix  and 
this  is  assumed  to  be  constant  for  each  vortex  and  of 
magnitude  determined  by  moderately  loaded  lifting  line 
propeller  theory(lfc).  The  equations  for  the  velocity  induced 
at  a  point  x,r,e  by  a  helical  vortex  of  unit  strength 
starting  at  }  ,P,  ,  are:(l&) 

oo 

ut=_L-  f  {r-j> cos  If'-Q))  * (x-jT;  s>sin  C 

4  f  [CK "  I')*"  -l-  r2- 1 -j>x  -lrj>  cos  f'J  3//z- 

oo 

U,  =  X-  f 

v  co$r']fe 


Live-0)  ^ostn.  0c-fO  COS  C ie' 
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The  basic  difference  between  the  propeller  calculations 
and  those  for  the  foil  is  that  these  equations  must  be 
integrated  to  obtain  the  induced  velocities  at  each  vortex 
due  to  all  the  others  rather  than  the  compar itively  simple 
equations  for  the  velocity  due  to  a  straight  line  vortex. 
This  Is  more  a  difficulty  in  computer  time  cost  than 
mathematics.  The  method  used  was  a  five  point  Gauss 
integration  performed  over  five  intervals  from  0  to  4380 
degrees  downs  t  re  am.  (16) .  The  integrations  are  performed  as  if 
the  trailer  began  at  the  downstream  plane  in  question  and 
extended  from  this  plane  to  infinity  downstream.  This 
simplifies  the  calculations  considerably  since  all 
integrations  are  performed  over  the  same  interval.  The 
results  are  then  multiplied  by  the  factor  analogous  to  that 
for  the  strip  theory  foil  model  on  page  14  of  Chapter  I: 

i+.  — ^ - 

1  fFTIT? 

to  retrieve  the  effect  of  the  portion  of  the  cause  vortex 
extending  from  the  plane  being  considered  back  upstream  to 
the  propeller  plane.  The  geometry  may  perhaps  best  be 
visualized  by  reference  to  Figures  (5.4), (5.5),  and  (5.6). 
The  velocity  induced  by  all  the  other  blades  on  the  sheet  in 
question  is  calculated  at  the  beginning  of  the  program  and 
assumed  constant  thereafter.  A  flow  diagram  for  the 
calculation  is  shown  in  Figure  (5.7). 

The  integration  of  the  velocities  at  point  X  to  obtain 
the  new  vortex  positions  at  plane  X  +A  X  leads  to  difficulty 
in  the  propeller  problem.  The  tangential  velocities  in  a 
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propeller  wake  (as  seen  from  a  blade)  are  often  four  times 
the  axial  velocity.  The  integration  of  the  Y,Z  components  of 
this  tangential  velocity  can  lead  to  a  spiraling  divergence 
from  the  correct  radius.  This  effect  is  corrected  by 
introducing  the  second  order  term  in  the  numerical 
integration  as  a  radial  velocity  as  follows: 

tr  ts  Lkz  error  maude  i-n.  sec. 

asin^  cl  -first  order 
n-umerlccLl  ilrit  ejection. . 
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Using  this  method,  a  calculation  was  made  on  a  20  foot 

diameter  propeller  with  a  elliptical  loading  along  the  span. 

* 

The  results  of  this  calculation  are  shown  in  Fig.  (5.1). 
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The  loading  for  this  propeller  was  extreme  <G=.lsin©)  and 
the  rollup  is  dramatic.  The  induced  velocities  for  this 
propeller  were  calculated  on  the  blade  by  calculating  the 
velocity  induced  by  each  element  of  each  vortex  of  the  sheet 
being  calculated  on  control  points  on  each  blade.  The 

velocities  induced  by  this  sheet  at  a  given  radius  on  each 
blade  are  then  summed  to  find  the  effect  of  all  the  sheets 
on  a  single  blade,  in  Figures  (5.2)  and(5.3)  the  velocities 
induced  by  a  blade  with  reasonable  loading  and  those  due  to 
this  extreme  loading  are  plotted  for  comparison  with  the 

results  from  moderately  loaded  propeller  theory  as  developed 
by  Kerwin(l6).  It  is  observed  that  again  this  method  gives  a 
heavier  loading  at  the  tip  than  expected  due  to  the  fact 
that  the  downwash  is  negative  at  the  tip#  increasing  the 

angle  of  attack.  The  results  at  the  hub  are  doubtful  since 
the  effect  of  the  hub  itself  is  neglected.  Work  is  in 

progress  on  a  scheme  to  include  the  hub  as  a  line  sink. 
Since  the  propeller  used  as  an  example  here  is  a  purely 
hypothetical  one,  a  purely  hypothetical  boundary  layer 

thickness  of  .05  feet  was  assumed. 

A  check  was  made  to  compare  the  contraction  of  the  tip 

vortex  trajectory  calculated  by  this  method  to  the 

experimental  results  of  Abraham  (18).  The  results  of  this 

comparison  are  shown  in  Figure  (5.8).  Abraham  photographed 

the  tip  vortex  (which  was  visible  due  to  tip  vortex 

cavitation)  and  plotted  its  radial  position  as  a  function  of 

c 

distance  downstream.  His  calculations  for  wake  contraction. 
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also  shown  in  Figure  (5.8),  did  not  agree  well  with 
experimental  results  but  the  present  calculations  agree 
quite  well  with  the  experiment.  The  inclusion  of  rollup  of 
the  trailing  vortex  sheet  is  apparently  important  in  making 
this  calculation. 


PACrfi  46 


X  induced  velocity  due  to  rolled  up  sheets 


FIG.  5.2  Acomparison  of  non  dimensional  axial  and 

tangential  velocities  for  a  lightly  loaded 
propeller  using  moderately  loaded  propeller 
theory  and  rolled  up  trailing  vortex  sheets. 
G  =  .03  sin©  ,  Hub  radius=».2R 


I  NOTATION  I 

(x,r,$)  INDl 
OF  ATRAIL 
TRAILER 

FIGURE  5 

: 


I 

! 
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FiG.5.7  LOGIC  FLOW  CHART  FOR  THE  LIFTING  LINE  PROPELLER  CALCULATIONS 
DESCRIBED  IN  CHAPTER  5. 
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CHARACTERISTICS  OF  MODEL  PROPELLER  TESTED  BY  ABRAHAM 
SEE  REF.  (18) 

Hub  radius  .1  ft.  Water  speed  13.35  ft. /sec. 

Propeller  radius. 5  ft.  Revs. /sec.  16.68 

TV  .301  4  Blades 

A  comparison  of  the  tip  vortex  radial  contraction  observed 
by  Abraham,  his  calculated  trajectory  using  Kerwin's  field 
point  velocity  program,  and  results  from  the  present  program 
appear  below: 


FIG.  5.8  A  COMPARISON  OF  RESULTS  FOR  THE  CONTRACTION  OF  A  TIP 


VORTEX  SHED  FROM  A  MODEL  PROPELLER 


CHAPTER  VI 


CONCLUSIONS 

It  has  been  shown  that  the  trailing  vortex  sheet 
problem  may  be  solved  numerically  for  an  elliptically  loaded 
foil  and  the  singularity  at  the  end  of  the  vortex  sheet  In 
the  potential  flow  problem  may  be  dealt  with  by  resorting  to 
an  inner  solution  employing  boundary  layer  theory  at  the  tip 
of  the  foil.  The  model,,  although  crude,  predicts  both  the 
minimum  pressure  in  the  tip  vortex  core  and  the  motions  of 
the  sheet  with  reasonable  accuracy.  The  use  of  a  strip 
theory  model  shortens  the  computation  time  to  the  point 
where  the  method  may  be  of  practical  use. 

The  lifting  line  propeller  problem  is  solved  for  the 
downwash  on  the  lifting  line  and  the  kinematics  of  the 
trailing  vortex  sheet  as  a  direct  application  of  the  foil 
solution,  a  rotation  of  the  coordinate  system  and  the 
inclusion  of  the  effect  of  other  sheets  from  the  rest  of  the 
blades  being  the  only  modifications  necessary. 

It  is  indicated  that  the  velocities  induced  by  the 
trailing  vortex  system  on  the  blades,  particularly  at  the 
tip,  may  be  quite  different  from  those  predicted  by 
moderately  loaded  propeller  theory.  This  may  lead  to  a  load 
distribution  on  the  blades  which  differs  considerably  from 
that  predicted.  Work  is  in  progress  on  the  inclusion  of  the 
lifting  line  model  in  the  lifting  surface  propeller  design 
method  described  by  Kerwin  (17), 
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APPENDIX 

COMPUTER  PROGRAMS 


imopucmN 

The  programs  used  In  this  work  are  listed  in  this 
appendix.  They  are  written  In  the  FORTRAN  IV  language  with 
the  exception  of  Professor  Milgram's  boundary  layer  programs 
which  are  written  in  the  MAD  language  for  the  M.l.T.  Time 
Sharing  System,  Flow  charts  are  included  for  those  programs 
which  were  not  flow  charted  in  the  main  text. 
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METHOP  1  COMPUTER  PRflQEAM 

This  program  is  the  iterative  solution  for  the  position 
of  the  trailing  vortex  sheet.  A  flow  chart  is  shown  in 
Figure  (Al).  This  program  uses  a  tremendous  amount  of 
computer  time  unnecessarily  and  it  is  hoped  that  it  will 
never  be  run  again.  A  listing  of  the  source  program  is  shown 
in  Figure  (A2).  Sample  results  are  shown  in  Figure  (1.2)  of 
the  main  text. 


ELLIPTICAL  LOAD  DISTRIBUTION  IS  ASSUMED 
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METHODS  2,3,4  COMPUTER  PROGRAMS 
These  programs  for  strip  theory  models  of  the  trailing 
vortex  sheet  from  an  elliptically  loaded  lifting  line  differ 
from  each  other  only  in  minor  details  as  indicated  by  the 
flow  chart  in  Figure  (A3).  A  listing  of  the  Method  4  program 
is  shown  in  Figure  (A4).  Sample  results  are  given  in  Figure 
(A5 ) . 


FIG. A3  LOGIC  FLOW  GRAPH  FOR  METHODS  2,3,4  OF  CHAPTER  I.  A  STRAIGHT 
LIFTING  LINE  WITH  ELLIPTICAL  LOAD  DISTRIBUTION  IS  ASSUMED 
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FINAL  COMPUTER  PROGRAM  FOR  A  LIFTING  LINE  FOIL 
This  program  includes  a  sink  at  the  tip  vortex  core 
center  and  calculates  minimum  pressures  as  well  as  motions 
of  the  sheet.  At  each  interval  downstream  the  velocity 
induced  by  each  vortex  element  on  selected  control  points  on 
the  blade  is  calculated  for  eventual  inclusion  in  a  lifting 
surface  method.  A  flow  chart  was  given  in  the  main  text 
(Fig.  3.1).  A  listing  of  the  source  program  is  given  in 
Figure  (A6)  and  sample  output  in  (A7). 
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FIG.A5  SAMPLE  OUTPUT  FROM 
METHOD  4 


FIG.A7  SAMPLE  OUTPUT  FROM 
AN  ELLIPTICAL  LOADING 
USING  THE  LINE  SINK  OF 
CHAPTER  3 
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LIFTING  LINE  PROPELLER  PROGRAM 
This  is  the  final  program  for  propeller  design 
purposes.  An  input  form  is  shown  in  Figure  (A8)  to  simplify 
use.  The  Fourier  series  coefficient  of  the  load  distribution 


are  the  coefficients  in  the  following  series  for  the 
circulation  on  the  lifting  line. 

Gr(Q)  -  K  (fl)  t  2.  Bn  s  in  (n.©)] 

n«i  J 

_  r  -  f^Hub 


SPAN 


where 


K=  irr  R  Uo 


SPAIN  r  a  -  Pvh.uk 

This  allows  the  application  of  this  method  to  other  than 
elliptical  load  distributions.  The  rest  of  the  input 
quantities  are  straightforward.  A  flow  chart  was  given  in 
the  main  text  (Fig.  5.7).  The  main  source  program  is  listed 
in  Figure  (A9)  and  the  two  subroutines  INVEL  and  INDUCT 
which  integrate  along  helices  for  the  effect  of  the  other 
trailing  vortex  sheets  and  the  other  vortices  on  the  sheet 
being  considered  respectively  are  shown  in  Figure  (A10).  The 
output  is  of  the  same  form  as  for  the  foil  problem. 
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Cma(I)-0.0 
OMIDH. 
otum-u. 

GlA(k)-)U. 

C.U(S)-S0. 

ttuCS)-120. 

GnA(7)-2kO. 

OU(U*)U, 

6lA(IWt«, 

QMAUQ)-k)20. 

C  CHANG E  TO  RADIANS 
OO  1000  H-1,10 

1000  G<U<N)-GiUE  4)-).lkl*t/ll0. 

OO  1021  nTHL-I.NOTrL 
1021  THENTKD-O.O 

n>uui» 

PITJ-PI-2. 

Pl-S.ikltl 

htwm 

CNOC  -PIIT-RPROP-VSHIP 

PluR-PI/$PAN 

NANO-NIUL*! 

C  MOHOCK  OF  VUKTICCS  IN  31 
0UM-NE3L 
01$P$l-Jt/©UM 

thick -otspse/2, 

C  0IST.  OF  FIR.T  TRAILER  In  FrOm  Ilf 
Ehqk-memq 

C  MO  WITH  FIN£  SPACING 

CNUL-k.-TMICK-ENON 
C  TUTAI  LENGTH  III  FINC  SPACING 
SPMF-SPAN-f  M>L 

C  TOTAL  LENGTH  WITH  MI0C  SPACING 
ttltKL-NOt* 

NSPS-NQTRL-l 

C  TOTAL  My  OF  SPACE  PjlHTs 

KSPSW-MUTRL-T-NCHU 
C  No  OF  SPACES  III  WIDE  SPACING 
0  T$  PSO  •  S  PNW/ F 10  AT  (  NS  PSJ> 

C  SPACING  M  WIOE  SECTION 

VTlL<l>-THICK»Ri*JS 
VUIEMuTKI)-SPAM-THICK-KiU.I 
C  FlhsT  TRAILER.  lii  FhiAI  CUO 

00  211  N-2.MCN0 
TTLL(N)-YIILEk-1)*UYSPsC 

211  VTLL<NOTkL>N*1)>SPA4-TTLL(N)  -RHU3  -rMU» 

C  OTHER  Thai LEh  POSITIONS  I  4  N.kHo.1  SPACI.4G 

t.iioll-span/2.  -rHUII 
YTII(EMuTRL"1)/2-1)-Yi4IULL 
c  Thai lek  at  hiosPa.4 

YSPS(l>-0.  -KHUJ 

C  SPACE  PUjlTlol  AT  CM  OF  SPA.4 

VS  PS  <  NS  PS) -SPA  II  -KHUN 

00  2  N-l.KCX) 

YSPS(N»2)-Y$P$(N)-0YSP*C 

2  ysp$ensp>-n)-span-ysp$Cn*i)  -khos  4Riiuo 

C  SPACE  positions  in  fi  <t  spacing 

Y  TlL<  MENU- l)-YTLl<NE -101 -Till  kK‘OYSPSW/2. 

C  FlA^T  TRAILER  PUSITIJi  II  UloE  SPACING 

VTllCNUTKL-KtND)-3PAN-YTllE.4C.40»l)  ‘KlIJS  -RHOO 
NTHAF-(WTKl-l)/2 

C  TKAILCM  JUST  TO  LEFT  OF  MIU.PA4 

nfakc-nend-2 

C  SECOND  TkAILEK  II  WiuE  SPAkC  SEC0I4)  $PAiC  POINT  I  I  HIUE  SPACING 
Do  212  M-nFAKE.NTHAF 
YTll<N)-YTLl(N-U*JYSPSI 
C  TRAILER  Position  1,4  WIDE  SPACING 

212  YTUCNOTKL-N.U-SPAN-YTIIEm)  -KMUJ  »RHUi 

NF  2 -l<$  PS/2 

C  SPACE  POSITION  TO  LEFT  Of  CE  HE 4 
00  21)  Hi, FAKE. NF 2 

c  space  positions  ii  wioe  spacing 

Y>PS(N)*YSPS(N-1)*UYSPS.I 

21)  YSPS(N$PS-N»2)-sPAr4<YSPS(N)  -KIWI  -RHU3 

DO  1  N- 1, MIPS 

C  GET  CIACULATIOI  AT  SPACE  POINTS 

1  GAHLI(I4)-CMC(YSPS(N>  -KHU3) 

GAHIL-CIRCOL  > 

CAHM-uAHOL/OuH 
Oo  %  M-l.MOTKL 

IFEVTU(N)  .LT.ESL  «KH03  )  >  Go  To  IS) 

CAjITkL(n)  — |.-IXIlL(i4»l)-GA.llL(,4>> 

GO  TO  % 

»)  oamtkieni-gamtk 

CAHIKLEN)  —  GAMTk 
k  Co  IT  I  NOE 

GCkUP-O. 

00  20)0  NFE-I.NIJL 
CCRUP-GCRUP* GA1TRLE  40TlL-l'NfE) 

1010  C04TINUC 

00  10)1  NFC-1.141  JL 

10)1  GAHTrLENUTKIM-NFCXCKUP/FIuATEMIRI) 

WRE  TEES,  100)  (GJkilKlOO.N-l.MiUO 
100  f  ORI-ATE  ///,  IX. 1  uAHTKL  */E  SF  10.S ) ) 

00  S  NTKl-l.MOTVL 
XEmTrL,  1)  -3.0 
YE  ITkL.D-TIILE  ITKI) 


MENTNL)-Y(nTKL.I) 

S  ZENTKb. l)-0.0 

NPOS-1 

C  WE  OUST  COT  THE  TRAILER,.  AT  ThC  LIFTING  LINE 
NCORE-E  Htl-n/2 
AC-Bl/2. 

C  COKE  RADIUS 

PI  2-PI  -PI 
A2-AC-AC 
Ak-A2-A2 
PITU-lt.-Pl 
PI  TI-I.pPI 
PIT2-I.-U 
PI2Tk-k.-PI2 
PI2T3-I. -PI J 
PI 2T1I-2. -PI 2T3 
CCKUP2 -CCHUP-GCKUP 

UU-Sl4KT<VSHIP»VSHIP*KPROP*KPRUP-U  lEGAvO,«EGA) 

Ul-W 

YU-0. 

PUP-PinF *Rikj»CCkup2/PI 2Tk  /A 2 
C  PkESSUkE  in  coke  AT  LIFTING  II «e 
FUP-PI  •A2-P|NF-)./Pim-KHO-GUCUP2 
C  FORCE  ON  UPSTKEAII  ENO  OF  COKE  AT  LIFTING  line 
S1-SPAN-3L 

nTuOUH-(IMTkl*1)/2-1 
C  HALF  THE  Thai  LEMS  -X 
WHfUlt.SlS)  PUP.FWP 

|)S  FOKMATE ///*  C0K|  PkESSUk-' FlS.k*  FORCE  On  UP$TKCh  ENO-'FIS.S) 

OO  US  NTKL-l.NOML 

US  M4I TEES. ISO)  NTKL.NPUS.XENT.4L.l),Y<NTfU,l).2(MTAL,l> 

C  FIRST  FI  NO  THE  STRENGTH  OF  THE  CAUSl  VORTICES  On  THE  OTHER  BLADES 
OIHCNSIUN  SPACE! 2 X).STK( 10). VCLAU*ES0).V£LTUIESO), VELKOBESU 
Oi)  2)kS  NTKL-l.MUTWL 
VELAO1(MTHL)-0.O 
VElTOB(NTML)-0.O 
VILKOIEnTRU-O.O 
2SkS  CONTINUE 

VSI NG-SPAN/10. 

RCAUSE 1)  -RHUB-VSI NG/2 . 

Oo  2kS0  H-2.10 

2kS0  RCaUS(H)  -KCAUS  (il-U*  VS  I  KG 
SPACEEU-O. 

OO  2kSl  W-2.il 

2kS 1  SPACEEH)-SPACEEM-1)«VSInC 
STkED— CIRCE  SPACE!  2) ) 

DO  2kS2  H-2.10 

2kS2  STkEH)— (CIRCE SPACEEN-1) ) -Cl kCEsPACCEM)) ) 

OO  2kS)  M-1,10 
CALL  INVELEW) 

VAilP-STREH)/Pt  Tk/RPROP 
OU  2kS)  NTKL-l.NOTNL 

VC  LAUt  t  NTKL  )  -VE  LAOBEnTRL  )*  VxMP-UAJB(rTRL.H) 

VEITOBENTKU-VELTOBE  NTRL)*VAHP*UTUiI(hTKL*M) 
VCW0IEhTRL)-V£IR03ENT«1)-VAHP-UR0BEnTRL,H) 

2k S)  CONTINUE 

c 

C  ROUTINE  TU  GET  INDUCED  VELOCITY  AT  BIASES 
C  LENGTH  OF  VUkTEa  CLEMENT  L 
REAL  U.LT.LZ 

01  HE  NS  I  ON  YKP!10.4).XKP(10.k),2KPa0.k).V*KP|10,S) 

2. VYKPE20.S). VZXPI  10.S)«NPLE  10).  VTKPEIQ.U 

c  No.acK  of  control  rui  its 
READES.IS20)  RPTuT 
NBTOT-KBTUT 
IS 70  FONMAT(U) 

WKITCES.U71)KPTUT 

SS71  Formate*  number  of  control  points  per  blaoi  *id 

REAO(S.tS72)  ENPl(O).J-l.KPTOT) 

»S72  format e xoi 2) 

WkI TEES. U2I)  (nPLE J). J-l.KPToT) 

IS7J  FORMATE  •CONTRUl  PUI'lTS  AT  SPACES  *  1012) 
c  find  control  points  on  first  blade 

OO  8S7k  0-1. RPTUT 
XKHEJ.l)  —  . 00001 
YKPEJ.D-YSPSEHPLEO)) 

VXKPEJ.D-O. 

YYKPEO,  1)  -0. 

VZKPE  J.  1)  *0. 
tS7k  COilTINot 

C  FIFO  CONTROL  POINTS  ON  OTHER  BLADES 
ALFN-0. 

ALFu-Pi  T2/F  LOATERBTOT) 

OO  UTS  NR-2.R0TOT 

AlFN-ALF  N-ALFU 

OO  UTS  RP-I, RPTUT 

XKPEKP.NJ)  — .00001 

FKPE  lP.NI)-TKP(KP,  DpCOSEALFN) 

ZKPERP.MlF)-YRP(  KP.U-SI  HALF  4) 

VXRP(KP,NB>-0. 

VYKPE  KP.NB  )  -0 , 

ViKPEKP,  N37-0. 

YTKPERP.NU-O. 

IS2S  CONTINUE 

00  ID  NPOS-2.NOPOS 
call  induct 

c 

C  thu  i *  the  big  Loop  for  each  uowm«tream  position 
WrI TEE*. 210) 

210  FoRMATESX.  1H  ) 

C  I  40UCCU  VELOCITIES  FRO, I  SIOT-SAVART 
DO  IS  NTNl-l.NoTRE 

C  For  each  control  vukTex 
C  INITIALIZE  VELOCITIES  AT  Till .  CP 
VELZENTKL.1)-0.U 
YlLY(NTKl,l)-O.U 
YELAENTkL.  d-o.u 
VELT  E  HTML.  1) -0.0 
V(LN(MlKL.l)-w.O 
DO  IS  HTkL-I.MOTRL 
C  For  each  cause  VORTEX 

IF  (IITRL-NlHL)  lk.  IS.  IS 
C  EXCEPT  ITSELF 

IS  CONTIMJC 

OY-YEHTKL.D-YMtKl.l) 

OZ-ZEMTRL. I)*ZEMTRI. 1) 

A-S(IRTE0T»DY-UZ-DZ) 

C  n  Is  For  the  CuUkui  poinTO 
C  M  IS  THE  VORTEX  IMUCX 

XBEQIN-O. 

xnu-xbecin-xentrl.i) 

Yr4*P-GA«4TRI(HTRI)/PI  Tfc/RPKUP-E  l.-XMU/ SORT  (A- A  -1NO-XNU)) 

C  IF  CP  IS  AT  THE  LAST  SHALL  SPACING  VORTEX 
IF  EnTRL-NEUI)  IIS.I1S.I1S 
US  IF EiITrE-EnENO-I ) )  I1S.U2.US 

•12  VA-VAHP-VAUVCnTKL.MTKL)-EUTSPSC*DYSP$N)/>./0YSP1W 
VT-VAlPtYruVENTRL.MTvi)-(OYsPSE-DYSP$R</2./0Y$PSN 
VH-VA)4P>VHUVINTRI,HTKL)-EUYSPSE*0YSPSN>/2./0YSPSW 
Go  To  111 
SIS  COHTIHOE 

If  CnTrL-ENCnUU))  t20.Sll.l20 
019  IFEHTRl-NCNU)  120.321.020 


FI&.  AS 


^ot  hephoducible 
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»I.**Hf.vru*d,T<L#HTMU«(UTX^l.u»SnJl/2./D»if$C 

vn.va»ip«v<uv(mTi<l.ktrumoT4P>c*..t..i>.,jj/2./o«s>'SI 

Cv)  U  tit 
<20  CuuTIwUl 

IHKT*t-(NuTXl-ll[i©)> 

12 1  It  (i>T41-(imTnl*(  h(NO-t)  >  Jtlk.Ul.ils 

*H  VA«VAaP*VAOY(NTHL.MT.il>*{uTSm«0fSi’SW>/2./OISPS( 
VUVAllP^TUYlNlRt.MlADMUYSRsEUTirSJJ/l./DYSPSt 
v<.»*4<««u*fiT*i.nTMV»«{ufsrse*0firtj)/2./t>^PSt 
Co  TO  III 
<2*  CuNTIMJt 

lf(wTNL'(ooT4L*(o(U0-t)))  121. <20. tit 
#20  If  |  NO  Ik  l  )  IJI.IJJ.42I 

DO  VA«VMf.YA0V<«*!»U.Nt.tD*{i;rsm.unfSH>/2./0T4»'SJ 
VT*VAlP*YTUV(*T>;L,MTRt>.CUY*Psl.U»SPS.i)/2./0YSVS- 
>TK.*fcl<«VK0Y<*.fKl..KTHlJ.tD»Sr5l.Ut4rtO)/2./0»;i»SJ 
U  Tu  111 

<21  VA*VAHP*VAOV(KTKl.HlHU 
VT*VAHP*VYu»(NT«l.HT».‘.> 

VR*YA*4P*V«UV<NTKl.Mt.*0 
111  TELAIkThL. t) *VCLA( XTmL, !)•»* 

«UTlNTKl(ll<mT(R?>a,|)4«! 

ULMNTm,  l)*VtU(NTNl,  1».*S 
is  Continue 

00  low  NTMl*l,MOT«t 

f  ACT04.1..XNU/S0-*t<l«'Ujf'f  10  VTU4TJT))..l.*r.y.XN0> 

VCU'NTNl.  l'«»tUO*Tfs,l).»Cl.A0|(«UU.f  ACTON.  VSHIY 
M»Tk0.i^»<rO,Tia.l,*YlUTKl..l).2<uTKV.X)»2<NTNl,l)> 

VE  LT(NTK«.„  M  *Vt LT (NtSt,  1 1  *  V[ L ToM N Tm L  ) •  4  ACTOk* M ( NT Al. ) 

Ull  V1U  l  NT*  Jl-YEUtNY  .I'-VELRJILNTAIWaCTuR 
c  mow  at/  yt.oc’m* 
e  mx*vti* 

oo  kTNL.i.MoJXL 

I^II/pEnTHU1  ^ 

1»U  m2<NVt,»>»VClt<NfRl.l>*VUt«l,t>/ft(H1rtl>*VClA<MTKl,U*z<NTKt.l 
c  on  Nt“  xi i  position* 

»  «*«,.  'l.woTia 

I'  ^nl,2i-IIVTRL.1>«01LX 
o/d  4Ux/r:>junTMi..i> 

"<  lTkL,2)*Tt  NTKL<l)*V(|lY(NTNL.  t)*wLL> 
''NTRl.2W<«Tl4L,l>*Vttl(l4YKl.l>*uttT 

4  l<NTRl>.Tl<lNTAL>*VtlT<NTAt.l>*UUT/R(*.T.IL> 

1  S«TlNUr;5,>  lt«l,}),T(  ItlUlMIuKl,!) 

ISO  fOAlATdlM.  20. S> 

C  NOW  THE  AXIAL  SuUCtJt 

C  INI  II  All  XI 

GCxOH-GCnUP 
GCRUN2-GORJP2 

C  STMT  WITH  UPSTREAM  CORE  STRENGTH 

VC*T(K>T><*1<HCU*C.|) 

IC*Z(N0T.4L*1-NCU.«C,1> 

00  210  NfCNI JL.NTRDUM 
X»-i(NJTKdl-Nf(,2) 

T*«TlH0Tn.|*tfl,2) 

02»A»S(IY-XC) 

OT.AJS(IY-TC) 

AfC(NFl)«SQKT(0Z*0Z*DT*9Y> 

**»UCI,1«I»  RFC(NF(),Nft 

IF (2V. Ll .XC. AMO. Y V.CI. YC)  THETAOFE  >.  O.U*ATAI4(OZ/OY> 

IFUV.IT.XC.ANO.YV.LT.YC)  THtTAlHft  >*  PI-  ATAN(UX/OY> 

If  (2Y.QT.2C.AN0.TY.tl.YO  iHCTAtNfC  >.*1.  AYA.4<UZ/»Y> 

If (XV.OT.2C.AiXJ.TV.CT.TO  T.UTAUFC  J-2..PI*  ATANJ02/0T) 

220  CONTINUE 

00  2X1  Nf(  •MAMG.NTliAX 
C  0(T  SOITENQCO  ANGUS 

AMC(MFt  >-(Tl«CTA(KX c  «l)*TilCTA(  .ft  •!> J/2. 

KnI T((4, 100DA <G(Hf (  J.Mfc 

221  C(OC(  Nf  {  )*£AHTNL(  -XJTRL*  l*Nf  C  >  •  AJi  (  A.<G(  Nf())/PIY2 
••NfC(MANC) 

C  riK!  VOATCX  ouTSloC  Coal  raj i os 
22S  P0-4l*Pli4f-CCK'UN2*-<,4U/PI2Tl/A2 
WA1TCU,110)>PON1 

C  MAST  cum  AT  PHtSsU«(,Ur.LV  CU*C  CIRCULATION 

PUUT*«.9 

POUT  AJUT*RJW*CCRUit2/PI  2T|*d  ./I/I-1./A7) 

0*11(1  »CCXON 

c  AIT  PRESSURE  out  TO  SHEET 
CO  221  MFC  ■  NAMG.MTHAF 
CAMIN*CA4tl!4*CtdC(Nf  (  ) 

222  POOT*POUT».tMO*OAHI.**]/PI2T**(l./Kft.(Nft*l>**2M./KfC(NFE>»*2> 

PUUT*PUUT-RMU*CA4II  J**2/PI  XTI/Rf  CL  NTH..F*  1  )<«i 

PONl*PONl»K)UT 

WMITCU.IIODPU  a 

C  PRESSURE  INCLUDING  OOTSlUt  SHEET  (Ff(CT 

I 0N1*PI  *A|*(  PI  Ml* POUT  ) 

»0*4i*f  JNl-RiU«CC-»D  12/PIT  U*Rjm«ClXO  42/PI  TI*A2*1  1  ./I/I-1./A2) 
NAITKi.lOOmUNl 
C  FORCE  ON  DOWNSTREAM  Cno 

Fl-foP*ONl 
WRITEU.lOIdd 
0JMy*I.»01-uu 

YN*AC/k./QllX*(-OUMV*$o4T(DU>IV*UOlV»k.«Fl/&i»o/P|/A2)) 

MWITKO,  10022VN 

C  IMOAKU  SI  XL  VELOCITY  IT  HO-lElfUM 
S1ML*P|T2«AC*VN 
WKIIEU.IOOOSI-XL 
U2*Ot*2.*VN/AC*OILX 
WNITEU.  II0SJU7 

C  AXIAL  VElgCllr  DUWNSTRE*!  XT  Cu.TI  XjlTY 
C  FIND  MOJ  HJOI  IS  POLUU  IN 
S2*S1»VN*0ELT 
Wn|I(U.lllO)S2 

Gsi*cuc<$n 

CSl-ClRCUX) 

MNtTEU.101IlGSl.Cs2 

OwCUK-GSX-GSl 

WAl  T(U,lll2>OXON.CAITrtL(40TNL*l«NANG> 

GCrD  4«<SCaUP«0oCuR 
WNiTEU.ll  lk)GCRO . 

CCaLM2*GCA0'4*i«CMD  i 
WNIT(U.101J)VO.V< 

I F  (aIUvj-yn)  -  .001  2)2),  22s.  220 

221  VU*VN 

C  TNT  A  NOTHIN  SI  XL  STRENGTH 
Co  To  72S 
22)  CONTINUE 

C  w(  HAVE  The  RIGHT  ST-EnuTh 
S1*S1 

OAHTNLE  wTnL*1*NA4G)«GAIInIL  wl  «l*  I  -  4A..w)  *0  .Co* 

»4**NdTRl* 1-NAmG 

WMI(U.1I1))CA.U«L(  <0VnL«!*II\<4).Mm 
CAITHU  40tNl»l-NCoNE)--AlTrtL(  .oT  *L*1-NConI )« j  .Con 
C  GET  VELOCITIES  Out  To  SI  a 

Do  220  Nf(  *l4L.4C.N(HAf 

VEIT(M|Ihl*  l*Nf  (  ,  J  )  *,  I  *.(Tv,.T<  »|T.  L»l*  «f  (.2>)/PlT2/NU(i4f()>.2 
2  21  V(U<WuUL«l-kf  £.  J  >*m»(  kiT-L*  l-Ml  ;.!)•(  «C*X<  40  Y  U*  l-ul  l.JI)/ 
l  (TC*t(  «JT.tL*l-Nf  (,2>> 


C  OCT  NEW  POSITIONS  INCLUUING  SINK  VELOCITIES 
00  222  HFt'NAMC.HTHAf 

1(.«OTNL*1''MF(,2)*T(NUTNL*1-NFC.2)*VLLT(i«JTNI*I>nF(,2)*OUT 
2(NOTNl*l-NFE#2>*X(NUTNL»l»Hf (,2)*VEIZ<NUTKL*I-Nf (,2>*0(LT 
WnITEU.IJOJNF  C.NPOS. X(>A>TKL* X'NF  (,  1 ) . T(  IOTRl*X-Nft,:},Z(MJTXl«l-N 
1FC.2) 

722  CUITIME 

WKITCU.D2)  VU.PUIl.fi 

1)2  FoaHATC  SI  «  VELUOTY"') IU.S*  PnESSU»E(LO/F TX-'f  10.S*  FORCE'FO.S) 
CCnuP-CCRDY 

C  INITIALIZE  FON  NEXT  P0*ITI04 
CCRUP2*0CKUP«CCNUP 
FuP*F0n1 
VW-O. 

IF(aFC(UAnO  .LT.  0)  NA4U-I1AI4G*) 

C  GET  VunTEx  CttHtwT  LENGTHS  FON  EFFECT  Oh  OLAO(S 
LX-OEIX 

0O  IS2S  Hol.MoTNL 
LV-T(M,2)-T(H,1> 

LZ*Z(.t,2)'Z(U,  1) 

c  get  s  FNun  control  points  to  vonTlx  eluients 
Oo  ISIS  NI-I.XOTOT 
00  ISIS  XP-l.KPTOT 
SX-UP(KP.NS)«(X(.l,l)/2.*X(.l,:>/2.) 
SY*TKP(XP,MO)-(V(M,l)/2.*V(H,2)/2.) 
si-jtfLir,Nx»-U(M,i)/i.»ic.i,:j/:.) 

S*SQNT(SX*SX*ST*ST*SZ*SZ) 

Coii«CA41TmL(H>/PI  Tl/S/S/S 

C  OCT  VELdClTT  INDUCED  IT  THI>  VunTCX  ON  THU  CUNTKL  POINT 
VXKP(KP,NI)*VXXP(XP,NI)*CUN*(LY*SZ*lZ*ST) 
VTXP(KP,NO)*VTCP(KP.MX)*CUN*(U*SX-LX*SZ> 
YXEP(KP,NI>*VXKPUP,N4)»lON*(U*SV-LY«SX> 

IS2S  CONTINUE 

00  221  NTKL'I.NOTXI 
XCHTML. 1)  •XLNTHL.2) 

T(NlRl,2)  •T(UTK1.2) 

Z(HTml.I)  •Z(NUL.Z) 

721  CONTINUE 
U1-U2 

WNITCU.72S)  02 

721  F04OUTUH  .’AXIAL  VELOCITY*'  FIO.S) 

)0  CONTINUE 

C  CALCULATE  TANGENTIAL  I  4DOCO  VELOCITIES 
DU  IS77  RP*l.RPTuT 

rkp>tkp(kp,ij 

OU  IS77  Nt-l.DTUT 

VTXP(KP,NI)*VZXP(KP,N3)*(TKP(XP,H«)/Hai>) 

1  *  VTXP(KP.I4t)*(  ZXPCKP,  ND/RAP) 

IS77  CONTINUE 

C  GET  NON  OIHENSIUNAL  AXIAL  Awo  TANGENTIAL  VEtUCIlT 
00  IS7I  XP*1,  KPTOT 
00  IS7I  NI*l.N|TOl 
VXKP(KP,  NI}*VXXP(XV,  NXJ/VCL 
VTXP(XP,Nll*VTILP(  KP.NS)/ VtL 
IS7I  CONTINUE 

C  GET  EFFECT  OF  ALL  ILAOES  u.  ONI  tLAUE 
00  IStl  XP*1. KPTOT 
VT-O. 

VA*0. 

DO  IStl  NJ.l.NtJlOT 
Vl*VT*VTKP( AP.NI) 

VA*VA*ViXP(  KP.NS) 

WAITCU.lSll )  TKP(KP,  D.VA.VT 
|SI2  FORMAT  (*  CUOlUS'FlO.S'AXlAL  NO  Vt  L’HO.S '  TANG' F  10.S) 

IStl  CUnTINje 

1001  fOM»AT(lSX'A4CLt»‘FlS  fc'AT  nTkL-'I). 

1002  fUNMATUSX'KUTAL-'FlS.S'AT  NTAL*'  I  )) 

100)  F  0K4LAK  XSX*  PUN1  FIRST  QUESS-'FIS.D 

100S  f  OMlAT(  ISX'f  UR.IAT  P0t41*'FlS.  *) 

109S  FORlATdSX'FUNl  I  41  Tl Al-’FlS.b) 

100k  FORMAT <XSX'Fl»'flS.k> 

1007  FOHMATdSX* VN*'f  IS.k) 

1001  FOR.IATdSX'51  4A*'F1S.») 

1001  FlM.IATd  SX'U2»'F1S.  *) 

1010  F0R.IATdSX'S2*'FlS.k) 

1011  F  OKMATdSX'GSl*'  FlS.  s'GSl^'F  IS.  k) 

1012  FORHATdSX'OGCJA-'f  IS.G’GAMTRLC  FlS.k) 

101)  FOXHATdSX'GMlTKLlNANGI-’f  lS.k'NANG»M  2) 

101k  F0R.UTdSX'GCKDN*'FlS.*> 

101S  fOR.l4Td)X'VU*'f  lS.k'V  4*'flS.k) 

CALL  EXIT 


F^G-.  (CO NT,} 
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SU3R0UTI  NE  IHVEL(.I) 

CU.IMUN  VAUV(  SO, 50) , VTUVC  SO, SO) ,  VRUV(50, 50) 

CUIVION  TH<  50 )  ,G<IA(  1 0)  ,RPitO P,GIAV(  SO  ) ,  PI  ,X(5O,2),Y(50,2),Z(50,2)  *  * 

Cu.vu.w  KUTUT, RCAUS(  10)  ,  UAOU( 50,  10),UTOB(50,  10)  , UK03 (50,  10 ) , NO  TKL 

DIMENSION  COEX(5),WEI  ullT(5),H(50)  ,COSKM(  1U ) ,  S  I.JKNC  10 )  i 

KT«K8ToT 

Do  1000  NTitL-l,NOTKL 
1000  R(  NTH L  )  -Y  ( flTIll,  1 ) /RPttU P 
DO  12  MTKL-1,M0TKL 
UAO8(NTRL,M)-0. 

UTOB(NTHL,H)-0. 

UROB(NTHL,M)-0. 

12  CONTINUE 
00  11  1-1,0 
OtLTA-GMA(  I  ♦  1 ) -G.  1A(  1  ) 

CU£X( 1 ) -. 0469 10-UELTA  t 

CoEX(2)-.23Q765*OELTA 

COEX( 3) -. 5-DELTA 

C0EX(4)-«  769235-DELTA 

CoEX( 5 ) ■• 95  3Q90-UELTA 

WE  I GHT(l)-. 059232-DELTA- 2. 

WE  1  GMT(  5  )  -WE  I  GUT  ( 1 ) 

WEI GHT ( 2) 11965  7 -DELTA- 2. 

WEIGHT(4)  -WEIGMT(2) 

WEI  GllT(  3  >  142222*Dt LTA*2 . 

00  11  1-1,5 

GNU-GtlA(  I  )-CO£X(  L)  , 

COSY-CQS(GUU)  } 

SI  IJY-SI  il(GNU) 

DO  11  K-2, KT 

CUSKN(K}-CO:>(2.*3.141S9/FLOATUT)*FLOAT(K«1)) 

SIHKN(K)-SlN(2.-3.141S9/FLOATUT)*FLOATU«l))  * 

cosinn-cosy-coskn(k)-si:iy-si.ik.ioc> 

SI  ill  iCN-S  I  NY*COSKN( K) -COSY-S 1  ilKN(K) 

VXI-GUV(.I  )«GI4U 
RZV-KCAUS(M)/RP«OP 
Oo  11  JCP-l,IWTia 

BOG-SURT((VXJ-VXI-K( JCP)--2  -kZV*KZV-2.-U(JCP)-KZV-C0SI KN) 

1--3  ) 

BOG-WEI GHT( L)/BUG 
EXGNU-VXI -RZV 

UAOU( JCP,M) -OAOB( JCP,M )♦ (AZV-RZV“RZV-K(JCP)*COSIKN)*BUG 
UT0B(  JCP,M)-UT0B(JCP,M)-(R(JCP>-GLRV(M)-RZV-GLRV(.O-C0SIKN-EXGNU*$ 
llNlKN)*  BOG 

UKOB(JCP,M)  -UROB(  JCP,M)-(GLi<V(M)-riZV*SlWI  KN-E  XGNU-CUSI  Kh)  *B0G 
11  CONTINUE 

return 

Eno 

SUBROUTINE  INDUCT 

COitMON  VAUV(  50,50),  VTUV(  50,50)  , VRUV(  50,50) 

CUUMOM  Tll( 50  ),GUA  ( 10)  ,RPKoP,GLKV( 50),PI,X(50,2),Y(SO,2),Z(50,2) 

COiVMON  KBTOT,RCAUS(  10)  ,UAOB(50,  10 ) ,UTJB( 50,  10 )  ,UKOB(50, 1Q),N0TRL 
OI.IENSION  CoEX(5),WEIGllT(5),R(50) 

Du  1000  NTRL-l,NOTKL 

1000  R(NTRL) -SQRT( Y  (  NTR  L,  1)-Y(  WTRL,  1 ) -ZOlTriL,  1) -Z ( NTRL,  1)  )/RPAUP 
00  11  NTKL-l,NOTRL 

DO  11  MTRL-i, NOTRE 
VAUV(NTRL,HTRE)-0. 

VTUV(NTRL,MTRL) -0. 

VRUV(NTRE,MTRE)-0. 

IF(NTKL-MTiU)  1001,  11,  1001 

1001  CONTINUE 
AHG-TH(I»TRE)-TH(MTKE) 

Oo  11  1-1,1 
0ELTA-GMA(  I ♦  1 )  — G< IA ( I  ) 

C0EX(1)-,  046910-DELTA 
CoEX( 2) 230765-OELTA 
CU£X( 3) -.5 •DELTA 
COEX(4)-.7G9235-OELTA 
CO£X(5)-. 95309  Q<*D£LTA 

WE  I GUT( 1)-, 059232-DELTA* 2. 

WCI  GHT(  5)  -WE  I  Gil T<  1  > 

WElGHT(2) 1 19657*OEL7A*2 . 

WE  I  GUT ( % ) -WE  I  GMT ( 2 ) 

WE I GHT( 3)-. 142222-DELTA- 2. 

DO  11  L-1,5 
GNU-GitA(  I  )*COEX<  L) 

COSY -CU$< GNU) 

S I NY-S I  (I  (GNU ) 

COSKN-COS(AflG) 

SI WKN-SI  N(  ANG) 

COS  I KN-CUSY-COSKN-S I NY-S I NKN 
SI  til  ilN-SI  NY-C03KN-C0SY-SI  NKN 
VXI  -G LR V ( . ITR L ) •  GNU 
RZV-U(ilT.CL) 

BOG-SURT (  (  VXI  -  VXI  -R( NTRL) -R(WTKL)-RZ V-kZV-2. -K( NTRL) -RZV-CUS I KN) 

1  —  3) 

BOG-WEI GHT( LJ/BOG 
EXGNU-VXI  -RZV 

VAUV( NTKL,MTK  L)  -VAUV(nTi<L,IITR  L  )♦(  AZV-RZV-RZV-R  ( NTRL) -COS  I RN) -BoG 
M-ltT«(L 

VTUV(NTKL,MTRL)-VTUV(NTRL JITt<L)-(K(NTUL)-GLRV(MT»a)-  RZV-GLiiV(M)* 

1  COS  I  KN-EXGNU-S I  Ml  .«N)  -BUG 

VKUV(NTRL,MTRL)-VRUV(NTRL,MT»tL)*(GL«V(MTRL)*RZV-SlNlKN-EXGNU-COSU 
1  II)  -BOG 
11  CONTI  NOE 

RETURN  * 

END 

R  19.210-6. 750 
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BOUNDARY  LAYER  PROGRAMS 

These  programs  are  listed  in  tne  MAD  language  purely 
for  the  sake  of  completeness.  They  were  written  by  Professor 
Jerome  Milgram  for  calculations  involving  seperated  flow  on 
yacht  sails.  See  Ref.  (11)  of  the  main  text.  The  listings 
appear  In  Figure  (All) 
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